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We consider a nonlocal version of the NJL model, based on a separable 
quark-quark interaction. The interaction is extended to include terms that 
bind vector and axial- vector mesons. The nonlocality means that no further 
regulator is required. Moreover the model is able to confine the quarks by 
generating a quark propagator without poles at real energies. Working in the 
ladder approximation, we calculate amplitudes in Euclidean space and discuss 
features of their continuation to Minkowski energies. Conserved currents are 
constructed and we demonstrate their consistency with various Ward identi- 
ties. Various meson masses are calculated, along with their strong and electro- 
magnetic decay amplitudes. We also calculate the electromagnetic form factor 
of the pion, as well as form factors associated with the processes 77* — ► 7r° 
and to — > 7r°7* . The results are found to lead to a satisfactory phenomenology 
and lend some dynamical support to the idea of vector-meson dominance. 



I. INTRODUCTION 

The Nambu-Jona-Lasinio (NJL) model |1]-J3| has long been used as a starting point 
for the description of light mesonic states as fermion-antifermion composites. It shares 
several features with QCD, notably a dynamically-broken chiral symmetry, with the pions 
as approximate Goldstone bosons. The model is based on fermionic fields interacting through 
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a local, chirally-invariant, four-point vertex. The local nature of the interaction produces 
a great simplification of the Schwinger-Dyson and Bethe-Salpeter equations. The main 
drawbacks of the model are, however, direct consequences of this locality. Specifically, they 
are that the loop integrals diverge (and so must somehow be regulated) and that the model 
is non-confining. 

Although the NJL model does contain regularization-independent information fUH and 
the results with various regularization schemes have been found to be qualitatively similar ||, 
the choice of any particular scheme lacks a sound physical motivation. A feature of many 
regularization schemes is that as well as the form of the cut-off, a particular momentum 
routing must be specified for loop diagrams with two-or-more quark lines [0. In practice 
a symmetric routing is often implicitly taken in order to maintain Ward identities. The 
regularization scheme must be specified yet further if one wishes to calculate beyond leading 
order in the l/N c expansion, a new cut-off being required for meson loops 

Another unsatisfactory aspect is that if low-energy theorems for anomalous processes 
(for example 7r° — > 77) are to hold then the anomalous diagrams must either be left unreg- 
ulated [[II]] or else additional terms must be added to the Lagrangian in order to respect the 
anomalous Ward identities |J . Related problems occur in the presence of vector interactions 
T(i| |13| if one attempts to apply the regularization prescription to both the anomalous and 



non-anomalous sectors. 



Many attempts have been made to generalize the original NJL model [14| to remove 
such unwanted features whilst retaining its successful phenomenological aspects One 
promising approach, which provides some motivation for the model studied here, is sug- 



gested by instanton-liquid studies fll5| . In that picture, the instantons induce an effective 



four-quark vertex, which is nonlocal but of a separable form. The separable nature of the 
interaction retains as far as possible the simplifying features of a local model, with the non- 
locality providing a cut-off on all loop integrals. A similar class of model has a separable 
dependence on the relative momentum of the qq pair and has been studied in Refs. ||TB|-|T7[ . 
Other models with simple interactions have been suggested based on various types of glu- 



onic field configurations within the QCD vacuum. For example, Efimov and coworkers fl8 
start with a constant (anti-) self-dual background field in Euclidean space and base their 
four-quark vertex on one-gluon exchange within such a background. Yet another recent 
model JB5J used a four-quark vertex mediated by a random colour-matrix, as an attempt to 
simulate a strongly-fluctuating background gluon field (see also pPfl ). We should mention 
that there are also studies of QCD Schwinger-Dyson equations based on one-gluon exchange 
forces between the quarks, often using effective gluon propagators |2l|j22l| (also see and 
references therein). 



Here we develop further the model proposed in Ref. f24|. This is based on a nonlocal 



separable, four-quark vertex. It is therefore similar to the instanton-liquid model of Ref. |T5 
the differences being that more general choices of the form factor and possible couplings are 
considered. The particular choice of form factor used in this model can lead to quark 
confinement, in the sense of a quark propagator without poles at real energies. It also 
ensures convergence of the quark-loop integrals, unlike that chosen in the separable model 



of Ref. [ . Only the pions and their scalar partner were studied in Ref. [ 23J . In the spirit of 



the extended NJL model 0-^JTyj2^-^D[ , we incorporate other mesonic degrees of freedom, 
such as the vector mesons. Including these particles enables us to probe the role of the 
confinement mechanism, since they have masses of around twice a typical constituent-quark 
mass. In a model without confinement, and with an otherwise reasonable constituent-quark 
mass of ~ 300 MeV, the p meson would lie above the qq threshold and could decay into free 
qq states. 

As well as the meson spectrum, we calculate hadronic decays of these particles and some 
electromagnetic processes, including p — > e + e~ and the pion form factor. With regard to 
the electromagnetic processes, we also examine whether there is any support for the concept 
of vector-meson dominance within a model of this type. 

In our present calculations, we work at leading-order in the 1/N C expansion. Our results 
are thus comparable to those of most NJL-model and Schwinger-Dyson approaches. At 
next-to-leading order, one has to include Fock terms arising from the four-point interaction 



and meson-loop contributions. Calculations at this order in the NJL model |||9| indicate 
that such corrections can be significant, especially in the scalar-isoscalar channel. This may 
explain some of the failures of the model to give a good description of amplitudes involving 
this channel. 

This paper is organized as follows. In Sec. |TJ we define the extended model and give the 
forms of the corresponding vector and axial currents. In Sec. |T| we present the forms of the 
quark and meson propagators and describe the means of coupling to external currents. The 
determination of the model parameters is described in Sec. |[V[ and results for hadronic and 
electromagnetic decays of mesons are presented in Sees. |V| and |V| respectively. 



II. THE MODEL 

Formally at least, one can imagine integrating out gluonic degrees of freedom to leave 
an effective action for QCD expressed in terms of quark fields only. As in the usual NJL 
model, we keep only two-body forces between the quarks, described by four-quark interaction 
vertices. Indeed, at leading order in 1/N C , all six-quark and higher interactions could be 
absorbed into effective couplings for the four-quark terms, by replacing extra ipYijj factors 
with their vacuum expectation values. This is just the procedure followed in the three- 
flavour extended NJL model with a six-quark, U( rebreaking 't Hooft determinant 



[ 31 1 . In the present work, however, we specialize to two-flavours with isospin symmetry. The 
action may be written as 

S = d 4 xip(x)(ifl - m c )ip(x) + / Y[d 4 x n H,^, x 2 , x 3 , x A ) ip(x 1 )T° : 'i/j(x 3 )'i/j(x2)Ti a ip(x i ). 

J i J n 

(1) 

Imposing SU(2)/jxSU(2) L xU(l)v symmetry requires that certain of the possible Dirac 
and isospin structures appear in the combinations: 

H X {1 ® 1 + z 75 r a ® z 75 r a ), H 2 (^T a ® 1(1 r a + 7 ^ 75 r a ® 7/ , 75 r a ), 

H 5 (r a <g> r a + t l5 ® z T5 ), tfeOv ® - a^ v r a ® a^r a ), (2) 



whilst the strengths of the following interactions are unconstrained by symmetry consider- 
ations: 



#3(7" ® 7m) > ^4(7 M 75 ® 7m7s)- (3) 

A variety of models of the above type can be found in the literature, differing in the 
ansatz chosen for {Hi(xi, X2, X3, X4)}. The original NJL model, for instance, has H\ ~ 
/ d 4 xl\ n S(x — x n ) and a constant coupling strength, whereas one-gluon exchange models 
take Hi ~ 8{x\ — x^)5{x2 — x±)D(x\ — x%). Our own approach is motivated in part by the 
instanton-liquid model Within the zero-mode approximation to that picture, there is 
a 2iVj-quark interaction, which is of separable form, 

Hi(pi,P2,P3,P4) = U 27r )' l Gif(Pl)f(P2)f(P3)f(P4)S(pi+P2-p3-P4)- (4) 



In the model of Ref. |15j the function f(p) is of a particular form and for two flavours of 



quarks the relation G\ = — G5 follows from the form of the 't Hooft determinant ||31|| . An 
interaction of tensor character is also present but is 1/N C suppressed. 



The model studied here is similar to that of Dyakonov and Petrov [ I5| , in that it is 
based on an interaction with the separable form (^J). However, we shall adopt a more 
phenomenological attitude towards the form-factor f(p) and the allowed couplings. Only 
interactions in the colour-singlet channels are considered in the present study. The G\ 
coupling (in the ladder approximation) produces the pions and their isoscalar scalar partner, 
a. The couplings in the spin-1 channels G2, G3 and G4 are responsible for the p, a 1; u and /1 
mesons. The G5 coupling allows the model to describe an isovector scalar and an isoscalar 
pseudoscalar meson. The lowest-lying meson with quantum numbers corresponding to the 
former is ao(980), whilst the latter is a non-strange state with the quantum numbers of the 
77 and 7/, which we shall refer to as 77*. 

We shall not consider the possible tensor interactions, described by the coupling Gq. As 
discussed in Ref. [J29|, these can contribute in the vector channels, where they give rise to 



anomalous magnetic-moment couplings of the vector mesons to constituent quarks. In the 



absence of any strong phenomenological need for such couplings we choose not to include 
them. 

For the sake of simplicity, all of the possible independent interactions (0,0) are assumed 
to contain the same form factor and to differ only in the constant coupling strengths, Gj. 
As in Ref . [JJ|] , we take the form factor to be Gaussian in Euclidean space: 

f(p E ) = exp(-p|/A 2 ). (5) 

This choice was shown to be able to give quark confinement. We have in fact examined 
the possibility of taking a different A for each of the independent couplings. Since this does 
not lead to any very significant effects, the results are not presented here. This is because 
the main features are dominated by the form of the quark self-energy which, in the ladder 
approximation, depends only on the G\ interaction. 

The usual, local expressions for the vector and axial currents do not satisfy the correct 
continuity equations when one uses the equations of motion derived from the action 
^). In order to obtain symmetry currents with the same divergences as the corresponding 
local currents in QCD, and hence to maintain the corresponding Ward identities, one has to 
introduce additional nonlocal terms in the currents. A Noether-like method of construction 



for these nonlocal pieces was developed in Ref. |[24|| . As a part of that procedure the following 
identity was used, 

5(x — xi) — 5(x — x 2 ) = / dX— — du5(x — z), (6) 

Jo dX 

z(X) being some arbitrary path from x% to x 2 - The divergence requirement determines the 
longitudinal component of a current which is, therefore, a path-independent object. In 



Ref. |M| the choice of path was irrelevant since the authors were interested only in the 
longitudinal component of the axial current, in order to determine the pion decay constant. 
The transverse part of the current, however, is sensitive to the particular choice of z(X). 
Indeed, ambiguity in the transverse current is a feature of any method used to construct 
a (partially) conserved current corresponding to a nonlocal action. When one wishes to 



consider electromagnetic processes, as in the present work, it is necessary to assume some 
form for the transverse current. This assumption forms an additional part of the specification 
of the model. In our present calculations, we adopt the straight line ansatz 



z(A) = (1 - \)x t + \x 2} (7) 

since this choice respects both Lorentz and translational invariance. In fact a number of 
the electromagnetic observables calculated in Sec.|V| are dominated by the local piece of the 
vector current and so should not be very sensitive to the choice of path. 

The nonlocal interactions in the model defined by (|l|-0) give rise to various nonlocal 
terms in the currents, whose momentum-space forms are presented below. (Note that where 
momentum derivatives with respect to Pi ±pj occur, then the combination ^fPj is under- 
stood to be held fixed.) For the isoscalar vector current, the nonlocal pieces all have the 
structure, 

J(!) = , 2 \ 12 Yl °i J II d ^Pn ^{pi)T^(p 3 )^(p 2 )Q ia ip(p 4 ) 5{pi+p 2 + q-Pz- Pa) 



i n 



x f d\f( P2 )f(pt) 9 f( pi + Xq )f(p 3 - q + \ q ), (8) 
Jo d{pi+p 3 )n 

which we refer to as type I. The sum over Gi(Tf (g)f^ a ) in (H) runs over the same combinations 
of couplings and Dirac and isospin matrices as those found in the action, (0||) . 

The isovector vector current also has nonlocal contributions of this type-I structure (^|). 
In this case the isospin and Dirac matrices appear in the combinations: 

G x (r a ® 1 + i 75 <g> i 75 r a ), G 2 (Y ® iuT a + 7^5 <g> 7^ a ), 

G 3 (Yr a ® 7,), G,(n 5 r a ® 7l /y 5 ), G B (1 ® r a + t l5 r a ® z 75 ). (9) 

Another type of nonlocal structure also arises in this current, 

Jfa = E °i J II d 'Pn ^(Pl)rtT b ^(p 3 )^(p 2 )Q ta T^( P4 ) 5( Pl +p 2 + q -p 3 -p 4 ) 



rl 

x I d\ 







d 

fipijfm) q, _ x /(P3 - q + M)f(p4 - M) 

d 

-/(Ps)/(P4) «7 ^f(Pi + Q - Ag)/(p 2 + Ag) 



(10) 



This type-II structure contributes in those interaction channels corresponding to isovector 
states. The Dirac matrices appear in the combinations: 

£1(275^75), G 2 (Y ® 7^ + 7^75® 7^75), G 5 (1®1). (11) 

Turning to the isovector axial current, the type-I terms are again present. These involve 
the matrix combinations: 

G ie abc (r c ® 2 75 r 6 ), G 2 ( 7 "75 ® 7.r a + 7" ® 7^), G^Ylsr" ® 7.), 
G 4 (Yr a ® 7,7s), G 5 e abc (z 75 r 6 ® r c ). (12) 

There are no type-II pieces in this current, but a third kind of nonlocal structure does occur, 

J (III) = / 2 Al2 H Gi / II ^Pn V'(Pl) r rV'(P3)^(P2)fi ia ^(p4) 5(P1 + P2 + 9 - P3 - P4) 



rl 

x / dX 
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/(P2)/(p 3 ) g , + x /(Pi + g - Ag)/(p 4 - Ag) 
.9 

-f(pi)fm)a? ; + Ag)/(p 3 - g + Xq) 

<9(p 2 +Ps) M 



(13) 



The relevant terms in this case are: 

d^r^l), G 2 e abc (Yl5T c ® lu r b ), G^®r a ). (14) 

It is straightforward to see that dependence on the path, A, does not appear in the 
longitudinal components of the currents. Since the form factor depends only on the square 
of its argument, one has in the case of type-I contributions 

d 1 d 

1 T^^ 1 + Ag)/(Ps - q + Xq) = -tt/(Pi + Ag)/(p 3 -q + Xq). (15) 

<9(pi+p 3 ) M 2 "A 

The A integral in q^J^ is then trivial, and produces a difference in form factors. Similar 
results can be seen to hold for the longitudinal components of the other nonlocal struc- 
tures (00). 



Useful checks on the above expressions for the currents are provided by various Ward 
identities which follow from (partial) current conservation [|3^]. In the case of the axial 
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current, an extension of the arguments in Ref. p4| can be used to show that the Gell-Mann- 



Oakes-Renner (GMOR) relation holds. For the vector currents, we have checked that 
the two-point correlator of vector currents is purely transverse, that the ^qq Ward identity 
is satisfied, that the pion charge is unity, and that the low-energy theorem for the anomalous 
decay n° — > 77 is satisfied. 



III. PROPAGATORS AND COUPLINGS TO CURRENTS 

A. Quark propagator 

The first step in the calculations is to construct the dressed quark propagator, by means 
of the corresponding Schwinger-Dyson equation (SDE). We treat it in the ladder approx- 
imation, which is equivalent to working at leading order in the 1/N C expansion. Fig. 1 
gives an illustration of the diagrams that are summed in this approximation. In terms of a 
momentum-dependent quark "mass" m{p) defined from the dressed quark propagator by 

S^ip) =j-m(p), (16) 

this equation can be written as[] 

roW ^ + iGl /^/||l±^|^ W . (17) 

The dressing at this order occurs only through the interaction in the isoscalar scalar chan- 
nel, G\. The integral here is very similar to that appearing in the quark condensate, apart 
from the extra interaction form factors. Just as in the original NJL model, the dynamical 
generation of quark masses is intimately connected with the appearance of a non-zero con- 
densate. In our numerical treatment of the model we evaluate loop integrals like that in ( |T7|) 
in Euclidean space, since the form factor has been defined for Euclidean momenta. Physical 



1 Here 'Tr' is used to denote a trace over flavour, colour and Dirac indices; the symbol 'tr' will 
later be used to indicate a trace over Dirac indices only. 
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results are then obtained by analytically continuing back to Minkowski space. Notice that 
the separable nature of the interaction produces a great simplification since the dependence 
on the external momentum p factorizes out of the loop integral. 
The solution to the SDE (|TTD can be written simply as 



m(p) = m c + (m(0) — m c ^jf 2 (p). (18) 

Hence to obtain the quark mass it is necessary to determine only the constant m(0). This 
can be done easily using iterative methods. In practice we choose to use equation (|17D to 
determine the parameter G\ for a given value of mo(0), the zero-momentum quark mass in 
the chiral limit. 

The denominator of the quark propagator, p 2 — m 2 (p 2 ), does not have a zero at positive 
(Minkowski) p 2 if m(0) is sufficiently large. This property provides a sufficient, although not 
strictly a necessary, condition for confinement. There are still poles in the quark propagator, 
but they are shifted into the complex p 2 plane. Such behaviour is not uncommon in models 
of quark confinement based on the solution of a Schwinger-Dyson equation P^J23| . pi 



in 



the ladder approximation. Because of the simplifications due to the separable interaction, 
the present model provides a convenient setting in which to investigate this mechanism for 
confinement. As pointed out by Lee and Wick |35j (see also [3f|), particles which have 



a complex mass of this type should not appear as asymptotic states if one is to have a 
unitary S-matrix. When amplitudes have been defined in Euclidean space, the prescription 
for analytically continuing them back to Minkowski space must respect this requirement, as 
described in more detail in Subsec. [Ill C| . 



B. Meson propagators 

The meson masses and vertex functions are found using the Bethe-Salpeter equation 
(BSE). For consistency with the SDE (0), we work in the ladder approximation. At this 
level, the separable interaction allows the qq scattering matrix, T, to be written in the form 
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T(pi,P2,P3,P4)=nf(Pn) S(Pi+P2-P3-PA)T{q), (19) 

n 

where we have denoted the total momentum of the qq pair by q = p\ — p 3 = p 4 — p 2 . The 
BSE, shown schematically in Fig. 2, may be conveniently expressed in terms of T as 

f(q) = G + GJ(q)f(q), (20) 

where G is simply a matrix of the coupling constants from the action (|1]-|4]) and J(q) is 
composed of the loop integrals 

Mi) = *Tr / ^f 2 ( P+ )f*(pJ)r i s F (pJ)r J s F ( P+ ), (21) 

where we have introduced p± = p ± \q. The quark propagators to be used in (^) are the 
dressed propagators obtained by solving the ladder SDE. 

The mesonic bound states are located at the poles of T. These can be determined from 
the equation 

det(l-GJ(qO) = 0. (22) 

The form of the action means that the matrix G is diagonal with respect to flavour and 
Lorentz structures. The scattering matrix, however, is only block-diagonal, since certain 
off-diagonal elements of J may be non-zero. In particular there is a loop integral that leads 
to mixing between the pseudoscalar and longitudinal axial channels. This nai (and rf f\) 
mixing is an example of the "partial Higgs mechanism" familiar from effective Lagrangians 
of 7r,p and a% mesons |37|j38[| . It leads to an axial as well as a pseudoscalar vertex function 



for the pion and rj*. With each flavour having an equal current quark mass, no such mixing 
occurs between the scalar and vector channels, as can be seen from the fact that the integrand 



in the corresponding element of fl2l|) is odd under p — ► —p. The absence of such a mixing 
means that the longitudinal vector channel is quite independent of the scalar one. It is 
therefore important to check that a pole does not develop in the former channel, since that 
would be unphysical. 
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For later ease of reference, we label the various elements of J as follows for the Dirac 
matrices inserted: 

J ss ■ 1 ® 1, Jvv ■ ® 7"), Jvv ■ <f 2 (-H ® 

^pp : «75®«75, </ap : m^{-%^ <g> 275), J Pj4 : m~ l (i^ b ® i$j 5 ), 

J T aa ■ T^(Yl5 ® 7^5), JjiA : m; 2 (-^ 75 <g> ^ 75 ), (23) 

where T MI , is the transverse projector, 

T, v = g, v - q -f. (24) 

In this basis, Jap = J pa- 

Near to the pole corresponding to a particular particle the relevant part of T may be 
represented as 

V(q) ® V(q) 



m 2 - q 2 



(25) 



where V(q) denotes the appropriate vertex function and any polarization indices have been 
suppressed. The vertex functions for the final and initial states are related by V = 7°V^7°. 
For the particles of interest, they are: 



V*(q) = (girqq - m~ 1 g 7rq J)i-f 5 T a , V a (q) = g aq 

Vps(q) = g P qqisT a , V ais (q) = ga iqq isl^ a , V ws (q) = g^qqis, 

V V *(Q) = - m^g^qj)^, V ao {q) = g aoqq r a . (26) 

For all particles except the pseudoscalars there is no mixing, and so each has a single coupling 
constant gi qq describing the coupling of an on-shell meson to quarks. These couplings are 
related to the corresponding loop integrals (|2"ID by 

Hqq 



gL dq 2 



(27) 



where S is the spin of the meson. The couplings of the pion to quarks, g nqq and g nqq , are 
given by: 
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, 1 - G 2 JUm 2 )) G 1 G 2 J PA {m*) 

y^qq U l D' (m 2 ) ' y^liy^qq D' (m 2 ^ ' 



where the prime indicates a derivative with respect to q 2 and the pseudoscalar-axial deter- 
minant D n (q 2 ) is defined to be: 

D^q 2 ) = (l " G x J PP (q 2 )) (l - G 2 J L AA {q 2 )) ~ G^J^q 2 ). (29) 

Similar expressions hold for the couplings of the rj*, with G 5 and G4 playing the roles of G\ 
and G 2 respectively. 

C. Loop integrals 

In Euclidean space, the loop integrals in the BSE (pl|) take the form 

Mq ) - -N c N f J (27r)4 {p 2 + + m 2 +){p 2 _ +m 2 _) ' (30) 

where is the appropriate Dirac trace and all momenta are to be understood as Euclidean. 
We have introduced m± here to denote the quark mass evaluated at p±. Consider this 
integral evaluated at some timelike momentum, q = (0,iq ). For a confining parameter set, 
each quark propagator considered as a function of energy has four poles at complex energies 
corresponding to a pair of complex-conjugate poles in p 2 . As go is increased these poles are 
translated parallel to the imaginary p 4 axis. For any given value of \p\, there is value of go for 
which the poles of the p- and p + quark propagators meet on the real p^ axis, pinching the 
contour of integration. For larger values of go the poles cross the axis and may contribute 
an imaginary part to the propagator in the meson channel, depending on the prescription 
used to continue the integral beyond the pinch point. 

The usual prescription based on Wick rotation of the integration contour would lead to 
an imaginary part of the meson propagator starting at the minimum value of go for which 
the contour is pinched. This would correspond to the opening of a threshold for decay of a 
meson into other states, and as described above (Subsec. [Ill A|) is inappropriate here. 
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One possible continuation, originally suggested by Cutkosky et al. [|36[| , amounts to the 



deformation of the integration contour displayed in Fig. 3. This prescription ensures that the 
resulting meson propagator does not develop an imaginary part above the pseudo-threshold 
where the contour is pinched by the complex poles, although it does mean that the prop- 
agator cannot be analytically continued past that point. Since the method is not unique, 
the choice of continuation prescription must be regarded as an additional assumption that 
forms part of the specification of any model with a quark propagator of this type. In our 



calculations we shall adopt the suggestion of Cutkosky et al., which was shown in Ref. |36 
to be consistent with the requirements of unitarity and macro causality. 



As discussed by both Cutkosky et al. |36| and Lee and Wick ||35|| , microcausality violations 



can occur in models with a Euclidean metric and states of complex mass. However, in 
order to measure such violations, Lee and Wick [[35| have estimated that one would need to 
create a wave packet of width <C 7 -1 , where the complex mass is M + ^ij. In any event, 
microcausality in the model is intrinsically broken by the use of an action with nonlocal 
interactions. 

In evaluating these integrals numerically, we take a contour in p 4 that runs along the 
real axis. For energies go above the pseudo-threshold, the prescription we use means that 
we have to add in contributions from the residues of the poles that have crossed the axis. 
Both the naive integral over Euclidean four-momentum in (|30|) and the residue contributions 
diverge at the pinch point, although these divergences cancel to leave a finite result |p6| . The 



cancellation occurs at the level of the integrated result rather than at all values of the three- 
momentum flowing round the loop. In our numerical work, we therefore need to regulate the 
two contributions when evaluating them separately. We subtract off a function which has 
the same divergence as the naive integral when p lies within a radius A of the pinch point. 
A similar function is used to cancel the divergent part of the residue contribution when \p\ 
is less than A from its pinch value and is chosen to cancel exactly with the piece that has 
been cut out of the naive integral. We have checked that our results are independent of the 
value of the regularizing parameter A. 
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It should be noted that the quark propagator of the present model has many complex- 
conjugate pairs of poles. Such an analytic structure is also found in the pion propagator of 



the N JL model within the proper-time regularization scheme [[JtJ • In our model these poles 
are present for both confining and non-confining parameter sets and their positions depend 
on the detailed behaviour of the form factor for large momenta. We therefore regard them 
as unphysical artifacts of the model. For the parameter sets considered here, the next set of 
poles would result in another pseudo-threshold at energies of ~ 2 GeV. Since the model is 
not intended to be credible at such momenta, these extra poles do not present a practical 
problem. 

D. Current couplings 

The electromagnetic or weak decay constant of a meson is given by the matrix element 
between the vacuum and that meson of the vector or axial current. In the present model 
one needs to include contributions from both the usual local current and the nonlocal pieces 
described in Sec. [TJ. The corresponding diagrams are shown in Fig. 4. 

We consider first the pion decay constant. The loop integral arising from the local part of 
the axial current is very similar to Jap-, except that only two (rather than four) form factors 
are present. One must also include a contribution from the Gi(i , ~f^T a ® 1) term with type-Ill 



structure (|T3|) in the axial current. As described in Ref. [p4 |, this piece can be written as 
a sum of terms each of which factorizes into a loop integral somewhat similar to that in 
the scalar quark condensate, together with a loop with a pseudoscalar insertion and a pion 
vertex function. It is convenient to refer to these loops as being one-quark or two-quark, 
according to the number of quark propagators involved. The contribution of the diagram to 

fn IS 

iGi f dH: Tr[^ + m(/c 2 )] r d 4 p TrV;(g)(jL + m_)7 5 r Q (^ + + m+) 

2ml J (2vr) 4 k 2 -m 2 {k 2 ) J (2tt) 4 (p\ - m 2 + )(pl - ml) AP+JAP-J 

x ~f\k)(f{ V+ ) + / 2 (p_)) - f{p + )f{p_)f(k)(f(k + q) + f(k - q))] . (31) 
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In the extended model there is another nonlocal contribution which is induced by the 
term Gii^f ® 7^75 T a ) with type-I structure (^j) in the axial current. The one-quark loop in 
this case has a vector insertion. Although the vacuum expectation value of i/j^ip vanishes 
by Lorentz invariance, a non-zero integral is produced by a combination of form factors that 
is anti-symmetric in the loop momentum. The contribution of this diagram to f w is 
-iG 2 f d 4 k TiYlfc + m{k 2 )} 



2ml J ( 2 4 k 2 -m 2 (k 2 ) 



[ f(k)(f{k + q)-f(k-q) 



v f d *P TrV;(g)(^_+m_) 7; ,7 5 r a (^ + + m + ) ^^ 2 
J (2^ (^-m 2 + )(p 2 _-m 2 _) 1 [P+)f [P - y {d2) 

These pieces of f w arising from the nonlocal current make significant numerical contributions 
and are needed in order to satisfy the Gell-Mann-Oakes-Renner relation [ 23]J52 |. 



A determination of the coupling strength of the a± particle to the transverse axial cur- 
rent requires the evaluation of very similar diagrams to those involved in f n . However, an 
important difference from the analogous nonlocal diagrams for /„. arises because the path- 
dependent integral over A for the transverse current is non-trivial. In this case a numerical 
integration over A is also required. With a type-I nonlocal structure (§), A is only present 
in form factors associated with one of the loops. This is not so for contributions induced 
by type-II (|10D or type-Ill (|T3| ) terms in the current, where the integrals for the one- and 



two-quark loops do not factorize. 

The couplings responsible for the electromagnetic decays of the vector mesons can be 
evaluated in a similar manner. Again the nonlocal contributions are numerically important 
and are essential if Ward identities such as 

^j^"{0\T{J^)O( 0m =0 (33) 

are to be satisfied. 



A common scheme to describe electromagnetic couplings in the literature f40j-P2|| , some- 
times called the impulse approximation, involves the jqq vertex only, neglecting irreducible 
couplings of the photon to more than two quarks. The ^qq vertex itself is chosen to be of the 



Ball-Chiu [^3j form. For the model studied here, use of the impulse approximation is not 
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an appropriate prescription for the calculation of electromagnetic observables. For example, 
as discussed in Sec. [VT|, it would not yield the correct value for the pion charge. In this 
model, the electromagnetic couplings are specified once a particular ansatz has been chosen 
for the nonlocal part of the vector current. The uncertainty inherent in the construction of 
the transverse part of the current is discussed in Sec. [TJ. Despite this, gross features of the 
nonlocal current would remain unchanged with different path ansatze. 

The dressed 'yqq vertex in this model is described in some detail below since it is a 
necessary ingredient in the calculation of many electromagnetic processes. An important 
constraint on the structure of this vertex, T M , is provided by the Ward-Takahashi identity 

^T ll (p,q) = S^ 1 (p + )-S^ 1 (p.), (34) 

where q is the photon momentum flowing away from the vertex and p is the momentum 
flowing through the vertex. (The isospin structure has been suppressed here.) 

The various pieces of the full 'yqq vertex within this approach are shown diagrammatically 
in Fig. 5. From the local current, there is simply a contribution to T M of the usual form, r y fJL . 
The nonlocal current induces contributions where there is a closed one-quark loop, similar to 
those in the pion decay constant discussed above. In the electromagnetic case, the diagram 
where the closed loop has a scalar insertion can be simplified by using the SDE fll7|) to 
express it in the form 

- (m(0) - m c ) [ d\^f(p + (A - \)q) . (35) 

Added to the local contribution, this would give the full vertex in a version of the model 
without vector mesons. Now, as a special case of Eq. (|i~5"D, we have that 

q ^ f ( p + (A " ^ = Tx f ( p + (A - (36) 

The A integral involved in gT M is therefore trivial, and so the Ward identity for the ^qq 
vertex ( |34"D can be seen to be satisfied by the sum of 7^ and (|35|). 

In the extended model with vector-meson degrees of freedom, there is another contribu- 
tion to involving a one-quark loop with a vector insertion, specifically, 
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-ilJ(p-)f(p + )G 2 N c N f 



d 4 k 



4k u 



Jo dX W f{k ~ q + Xq)f{k + Xq) - (37) 



(2vr) 4 P -m 2 (k 2 ) 

In addition there are pieces which contain the propagator of an intermediate qq state, as 
constructed from the ladder BSE (|20l), in the vector channels. The contribution of the 
longitudinal channel to is 



i-of(p-)f(p-\ 



G 2 N C N 



f 



d 4 k 



f(k-)f(k, 



l-G 2 Jb v (q)J (27r) 4 (k 2 _-m 2 _)(k 2 + -m 2 + ) 

xtrf M (A;, q)(]t- + m_)^(# + + m+), 



(38) 



while the transverse channel gives 



Q 



f(P-)f(P-\ 



G 2 N c N f 



G 2 Jyy 



(?) 



d 4 k 



f(k-)f(k, 



(2tt) 4 {k 2 _ - m 2 _){k 2 + - m\) 

xtrf q)(jt- + m-)j v (ft+ + m+), 



(39) 



where V J '(k,q) is the two-quark-irreducible ^qq vertex consisting of the sum of 7 M , 
and Q3"7|). In these expressions, m± denotes the quark mass evaluated at k±. 

To check that the additional contributions (P7|-|3"9"|) in the extended model are consistent 
with the Ward identity (0), we note first that the quark propagator is unchanged and 
so the sum of the local piece and ( |3~5|) still saturates the identity. In the contribution 
of expression (|37f) to gT^, Eq. ([15]) enables the integration over the path variable to be 
performed. This part of gT M is then 

d 4 k 4q ■ k 



-i\f{p-)f{p + )G 2 N c N f J 



-f(k)(f(k + q)-f(k-q 



(40) 



q* 'J (2vr) 4 k 2 -m 2 (k 2 ) 

The purely transverse piece (|39|) , which involves a propagating p-meson, obviously does not 



contribute to the Ward identity. Thus cancellation must occur between ( |40[) and the piece 
coming from ([38]). To demonstrate this explicitly, we make use of the fact that gT M is given 
by the sum of (|40|) and the right-hand-side of (|3j). The contribution to q^T^ involving the 
longitudinal qq intermediate state (|38|) can then be expressed as 



G 2 N c N f 



q 2 1 - G 2 J^ v {q) 



f(P-)f(j>-\ 



-G 2 J^ v (q) 



d 4 k tr(^ + m_ — m + )(f,_ + m_)^(^ + + m_| 
(2^)4 (4 - m ^)(F -m^) 

d 4 fc 4g • k 



-f(k+)f(k. 



(2n) 4 k 2 -m 2 (k 2 



-f(k)(f(k + q)-f(k-q) 



(41) 



The Dirac trace in the first line of the above expression may be written as 

A(q ■ k„)(k 2 + - m 2 + ) - A(q ■ k + )(k 2 _ - m 2 _). (42) 

Hence in each of the resulting terms, one of the factors k\ — m\ can be cancelled in the 
denominator of the integral. Shifting the integration variable from k to k± as appropriate, 



the first term inside the square brackets of fl4l|) may be rewritten in the same form as the 



second, demonstrating the required cancellation. Note that the discussion above has referred 
to the presence of the G2 coupling in the isovector interaction channel. The results in the 
isoscalar channel are completely analogous, with the replacement of G2 by G3. 

For the purpose of practical calculations, it is convenient to collect together the various 
contributions to the vertex in the form 

r„(p, q) = IvQ + (7, - f(P-)f(p + )B(q 2 ) 

-2Q £ d\(p + (A - \)q),m\p + (A - \)q), (43) 

where the prime denotes a derivative with respect to the square of the momentum argument 
and we have reinstated the flavour structure using the charge matrix Q = |(r 3 + |). The 
function B(q 2 ) accounts for the presence of vector mesons in the model and has the form 
B(q 2 ) = ^(T 3 B 2 (q 2 ) + |-B 3 (g 2 )) where the functions Bi(q 2 ) are 

(q-k) 2 \\ 



Am_m + - Ak 2 + q 2 + - [k 2 

3 V 1 

q 2 



~{m + + m_) (k 2 - j £ dX m '(k + (A - \)t 



(44) 



and the Ai(q 2 ) in the above equation originate from the one-quark loop with a vector inser- 
tion and are given by 

Jz_ f d A k k 2 - (q- k) 2 /q 2 



A(q 2 ) = -jGiN c N f J 



(2tt) 4 k 2 -m 2 (k 2 ) 
jT 1 dX(f'(k + Xq)f(k -q + Xq) + f(k + Xq)f'(k - q + Xq)). (45) 
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Writing Ai(q 2 ) and Bi(q 2 ) in Euclidean-space form and performing an integration by 
parts on one finds that -Bj(O) = 0. This is simply a consequence of the differential form 
of ©, 

r,(p,0)=Q^Ls F \p). (46) 

Hence, in evaluating processes where the photon is on-shell, the 'yqq vertex is unchanged by 
the existence of vector-meson degrees of freedom in the model. 

For the coupling of an on-shell vector meson to the vector current, one evaluates a quark 
loop between the meson vertex function and that part of the jqq vertex which does not 
include the contribution from the propagating transverse vector channel, fl3~9"|). Attempting 
to include such a diagram would merely amount to the addition of another bubble onto the 
vector-meson chain. Since (|39"D is purely transverse, the Ward identity for the vertex still 
holds. 



IV. NUMERICAL FITS 

The model as defined in Sec. [n] above has seven parameters: the current quark mass 
(m c ), the range of the form factor (A) and five possible coupling constants. Considering first 
the couplings G\ and G2 only, we have chosen to fit the quantities: m n = 140 MeV, f n = 93 
MeV and m p = 770 MeV, which do not depend on the remaining three couplings. This leaves 
one parameter undetermined which we have taken to be mo(0), the zero- momentum quark 
mass obtained in the chiral limit of the SDE (|I~7D. Each of the other three couplings may 
then be independently fixed to reproduce the mass of the corresponding meson: G3 is set by 
requiring = 783 MeV ; G4 by = 1282 MeV and G5 by m ao = 982 MeV. The meson 
masses are determined from Eq. fl22"|) , whilst f n is given by the coupling of the pion to the 
axial current. The contributions to f n from the nonlocal pieces of the current are significant: 
the scalar and vector loop terms described in the previous section are respectively ~ 35% 
and ~ —10% of the total value. 
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In terms of m (0), the possible fits have a restricted range. Having a coupling strong 
enough to realize confinement requires that m (0) > 270 MeV. Below that value, the model 
should only be used up to an energy corresponding to the appearance of the qq continuum 
at twice the value of the (purely real) quark mass. In fact only a very limited range of non- 
confining sets are possible because the empirical masses of the vector mesons are located in 
this continuum for m (0) < 250 MeV. 

An upper limit on the acceptable values for rao(0) is imposed by the behaviour of the 
meson propagators above the pseudo-threshold. The dramatic changes of behaviour which 
can occur beyond this point may be seen in Fig. 6, where the denominators of the propagators 
in various scattering channels are plotted for a parameter set with m (0) = 300 MeV. In 
this case the pseudo-threshold occurs at an energy of 895 MeV. For larger values of m (0) 
this energy decreases. As suggested by the behaviour in Fig. 6, for large enough m (0) two 
additional poles appear in the transverse-vector channel above the p pole. The first of these 
extra poles has a residue of the wrong sign to describe a physical particle. Although one 
might hope to consider parameter sets with the extra poles, provided that they lie well above 
the energies of interest, in practice this is possible only for values of too(O) within a very 
narrow range, ~ 320 — 330 MeV. 

A pronounced change in behaviour beyond the pseudo-threshold also occurs in the 
longitudinal- vector channel and seems to be important in ensuring that no poles are present 
in this channel. An unphysical pole does occur, however, in the pion propagator (Fig. 6). 
This unwanted pole is located between 1.3 and 1.45 GeV, depending on the parameter set 
used. One should therefore only attempt to use the model at smaller energies. Although a 
similar pole also appears in the 77* propagator, it lies at a higher energy than in the pion 
case. 

We present here numerical results for parameter sets which lie near each edge of the ac- 
ceptable range for m (0). The variation of results over the full range is generally monotonic; 
where it is not, the dependence on m (0) is fairly weak. Specifically, we have chosen to 
work with m (0) = 280 MeV (henceforth referred to as set A) and 320 MeV (set B). Details 
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of these parameter sets are given in Table 1. Values of the zero-momentum quark masses 
calculated for these parameters at non-zero m c are also quoted in Table 1. 

Including a non-zero m c in the SDE has a significant effect on the solution, a current 
mass of ~ 10 MeV leading to an increase of ~ 50 MeV in the zero-momentum dynamical 
quark mass. It is thus worth examining the deviations of /„. and from the values which 
would be obtained at leading order in the chiral expansion. If one sets m c to zero and keeps 
all other parameters fixed, then the corresponding values for the chiral limit of f n are 84.6 
MeV for set A and 85.1 MeV for set B. To leading order in m c , the pion mass is given by 
the GMOR relation 

m 2 Jl = -m e (fop) , (47) 

where f n and the quark condensate are evaluated in the chiral limit. This expression gives 
values for m n of 143.6 MeV for set A and 143.2 MeV for set B. Comparing these results with 
those of the full calculations for these parameter sets, we see that higher-order terms in the 
chiral expansion contribute about 10% to f n . In contrast, the corresponding contributions to 
m n are less than 3%. These observations suggest that the restoring forces against deviations 
from the chiral circle are rather weak in this model. One would therefore expect to find a 
light a meson in connection with this softness of the vacuum. This is indeed the case, as 
will be discussed shortly. 

In the chiral limit, the model quark condensate is — (206MeV) 3 and — (189MeV) 3 for 
sets A and B respectively. With non-zero current-quark mass, the condensate integral is 
quadratically divergent. If it is regulated by subtracting the perturbative condensate, slightly 
higher values of — (212MeV) 3 and — (193MeV) 3 are obtained. These are similar in size to 



values for the condensate determined from QCD sum rules jP|,fEJ. However one should 
bear in mind the fact that the condensate in QCD is a renormalized and scale-dependent 
quantity and so one ought to be careful about comparing it directly with the value obtained 
in a model of this type. 

Table 2 lists the positions of the first few sets of poles in the quark propagator. Since we 
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only consider the first group of poles to have physical relevance, the model should be used 
only for energies of less than twice the real part of the second set of poles. This limit is 2.3 
GeV and 1.9 GeV for the parameter sets A and B, and so is sufficiently far above the upper 
limit imposed by the unphysical pole in the pion channel not to be of practical concern. 

In Table 3, the calculated meson masses are given, along with their on-shell couplings to 
quarks, as defined in Eqs. (|2~7h f29"l) . The scalar isoscalar state is rather light. For comparison, 
the mass of the corresponding particle in the NJL model H,[|^7| is m 2 = m 2 + 4m 2 , where 
m is the mass of the constituent quark. Interestingly, the a mass in our model varies only 
slowly with the dynamical quark mass. 

The empirical determination of the mass of the corresponding scalar isoscalar meson 
is still unclear because of the very strong coupling between that meson and the two-pion 



channel. While some analyses find masses of the order of 1 GeV [46], others result in much 



smaller values fl"7 |. The a masses obtained in this model, like those in the NJL model, 



are compatible with the latter. However, before taking them too seriously, one should 
remember that the strong coupling to two pions also means that 1/N C corrections are likely 
to be important in this channel [§,0. It is thus premature to draw conclusions about the 
ability of the model to describe this channel until such corrections have been calculated. 



The calculated a\ mass is somewhat smaller than the observed 1230 MeV [jiql . In the case 
of parameter set A, it lies a little below the pseudo-threshold, but for most of the range of 
admissible m (0) it is above that energy. The p — a\ mass splitting increases with increasing 



constituent quark mass, although not so rapidly as suggested by the NJL |27| expression, 
m ai = m p + 6m 2 , obtained from the derivative expansion of the bosonized model. As a 
consequence of the upper bound on the constituent mass, which follows from the effect of 
the pseudo-threshold on the transverse-vector channel, it is not possible to simultaneously 
reproduce the empirical values for both the p and a\ masses in the ladder approximation. 

Since there are important flavour-mixing effects in the isoscalar pseudoscalar sector, a 
realistic calculation for these mesons would require a three-flavour version of the model. 
The if mass in the present two-flavour model should not therefore be directly compared 

23 



with experiment. It is nevertheless reassuring that this mass lies between the physical 77 
and 1]' masses of 547 and 958 MeV respectively. Another possibly important feature in the 
description of the state is the axial-pseudoscalar mixing with the longitudinal f\ channel. 



Indeed, in a Bethe-Salpeter study of a three- flavour model [EJ, the ^75 term in the vertex 
function was found to make significant contributions to the mass (~ 70 MeV) and decay 
constant of the 77. A similar effect has also been found in the NJL model |[29|| . In the 
present model, if the f\ is omitted by setting G4 to zero, then the 77* mass for parameter 
set A is reduced by around 20 MeV, whereas for set B it falls by over 110 MeV. These 
rather different behaviours are another consequence of the dramatic changes in the meson 
propagators which can occur at the pseudo-threshold. When G 4 = 0, the 77* mass lies below 
the pseudo-threshold for the full range of parameter sets. For non-zero G4, the mixing 
increases the 77* mass and for parameter sets with m (0) > 310 MeV the mass is pushed 
above the pseudo-threshold, where the effect can be greatly enhanced. In addition, for these 
parameter sets the gradient of the determinant D v * (cf. Eq. (|29"D) changes significantly above 
the pseudo-threshold and so the coupling of the 77* to quarks is considerably stronger. 



V. HADRONIC DECAYS 

At leading order in 1/N C , the three-meson vertices are calculated from a quark loop 
with insertions of three vertex functions. We have obtained several physical decay ampli- 
tudes from such diagrams with all of the mesons on-shell, where the vertex functions are 
unambiguous. 

For an initial state of momentum q decaying to particles with momenta q% and q2, the 
quark propagators in the triangular loop are evaluated at p ± \q and p + \{q% — q2). If 
the initial state has a mass which is greater than twice the real part of the quark pole, 
then its decay modes will be sensitive to pseudo-threshold effects. By analogy with the loop 
integral in the BSE for that particle, residue contributions must be taken into account in the 
three-point diagrams. It is also possible that further residue contributions may be required 
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if a final-state particle lies above the pseudo-threshold energy, but such a situation does not 
occur in practice for any of the amplitudes considered. 

The couplings that we have considered are defined below, and their numerical values are 
given in Table 4, along with the corresponding decay widths. 

(7r a ( gi )7r 6 ( g2 )|<7(g)> = -g^J ab 
(^(<?iK(<72)|p a (g)> = ig P ^e ab \q 2 • 6 - Ql ■ e) 
(a(q 1 )7r b (q 2 )\al(q)) = ^ig ai aJ ab {qi -e-q 2 -e) 

</A<7iKfe)K(g)> = e abc (g aipv (e; ■ e a J - h aip7T (q 2 ■ e* p ){q 2 ■ e a J) (48) 

If the sigma meson of the model is to be interpreted as corresponding to the scalar 
particle of the linear sigma model then its coupling to two pions should be strong. Whilst 
the values in Table 4 do not indicate a particularly broad state, this width is significantly 
reduced by the small available phase space. To provide a useful comparison, we consider 
the prediction for the two-pion coupling g U7Tn from the linear sigma model. In that model 
the coupling is g U7Tn = — m^)/ f n which, for the a masses of parameter sets A and B 
gives g CT7r7r = 1901 MeV and 2122 MeV respectively. These values are ~ 30% larger than 
the results for our model, indicating that the coupling to pions of the scalar meson in the 
nonlocal NJL model is qualitatively similar to that in the linear- sigma-mo del. As already 
mentioned, this strong coupling indicates the importance of going beyond leading order in 
1/N C . 

The p meson decay width compares reasonably well with the empirical value of 151 MeV. 
In contrast, the corresponding leading-order calculation in an extended NJL model, using 
the physical p mass, significantly underestimates the decay rate |5(|. We find that it is not 



possible to simultaneously reproduce the empirical mass and decay width of the p at leading 
order in our model. However, if the model parameters for a given mo(0) are refitted to the 
p width rather than to its mass, the results for observables are not qualitatively different 
from those of the original fit. For instance, this procedure would increase the p mass itself 
by ~ 20 to 60 MeV. 
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The coupling g ai<J7T is not a direct observable, although the process it describes would be 
involved in the physical decay of a\ — > 3it. The partial widths for a\ — > an found in this 



model are similar to those in other extended chiral models |5T|J52| . In contrast, the Particle 
Data Group |48[] quotes an experimental upper bound on the final state 7r(7T7r)s of ~ 0.7% 
of the total a\ width of ~ 400 MeV. The strong couplings obtained here suggest that the 
model may not be consistent with such a result. However, the situation is far from clear. 
The two-stage process a\ — > an — > 3n would have to be integrated over the momentum of 
the intermediate scalar meson. For momenta corresponding to a virtual a, the a\an and 
aim couplings may be reduced from their on-shell values. A hint that this may indeed be 
the case is provided by the ann loop integral, which vanishes at a total energy of around 
800 MeV. Furthermore, there is also an amplitude for the (7T7r)s7r final state originating 
from a direct, four-point ai — > 3n diagram which we have not yet calculated. Although an 
estimate of this coupling in the derivative expansion of the extended NJL model p8| finds 



it to be small, it is conceivable that it might conspire to cancel some of the amplitude due 
to the intermediate scalar meson. However, as discussed in Ref. |53 |, such a cancellation is 



not required by any underlying principle such as chiral symmetry. Finally one should also 
note that, as already discussed, leading-order calculations in l/N c may not be sufficient to 
describe a process in which the scalar isoscalar channel plays a significant role. 

The dominant decay mode of the a>i is to pn. Although the parameter set B does give a 
credible, broad width, with set A the state seems to be very narrow. Since the final state has 
a combined mass of 910 MeV, the allowed phase space for the decay is drastically reduced at 
the model a\ masses as compared to the empirical mass. Using parameter set A, the a± mass 
is only 946.8 MeV and so the small decay width of 44 MeV may simply be a consequence of 
the phase-space suppression. 

In order to examine whether the ciipn coupling is reasonably well described by the model, 
we compare it with the description of the same process using a phenomenological mesonic 
Lagrangian. The CCWZ formalism originally developed in Ref. |53J is particularly conve- 
nient for this purpose since the a\ mass can be set to any desired value without violating the 
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constraints of chiral symmetry. A suitable Lagrangian for the comparison is one obtained by 
converting the simplest Lagrangian in the massive Yang-Mills scheme [[J7[] into its CCWZ 
equivalent and then adjusting the a\ mass.Q Details of the relevant change of variables can 
be found in Ref. j38]: the appropriate Lagrangian is equation (5.6) of that paper. We take 



f n = 93 MeV, m p = 770 MeV and Z 2 = 1/2 as the parameters specifying the original 
massive Yang-Mills Lagrangian. For the empirical a\ mass of 1230 MeV, this gives a broad 
state of width 490 MeV. When we use the a\ masses found in the model with parameter 
sets A and B, the effective Lagrangian gives significantly smaller widths, 23 and 132 MeV 
respectively. This suggests that the small widths found in the model are largely due to the 
the small a x mass rather than any underestimate of the coupling strength. 

The amplitude for the decay ai — > pir is a mixture of s and <i-wave components. In terms 
of the decay parameters defined in (f48|), the ratio of d to s-wave amplitudes is 

(49) 



(E p - m p )g aipn + | 


%\ 




(E p + 2m p )g ai p n + 


K 





where E p and q are the energy and three-momentum of the p, in the a x rest frame. This 



quantity has been determined by the ARGUS collaboration |55j from r-decay data to be 
—0.11 ± 0.02. The effective Lagrangian approach discussed above requires higher-order 
couplings in order to obtain a non-zero h aifm and so this ratio provides a test of such higher- 
order effects. From the values for R in Table 4, we see that, for the parameter set A, the 
ratio is rather low, but the value for set B is consistent with the observed one. Overall, set 
B provides the better description of both the a\ mass and its hadronic decays. 

VI. ELECTROMAGNETIC PROCESSES 

Further tests of models where mesons are constructed as qq composites are provided by 
electromagnetic processes and form factors that probe the internal meson structure. For this 



2 A similar conclusion is reached if one starts with the simplest Lagrangian for ir, p and a± in the 
hidden-symmetry formalism. 
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it is necessary to specify the photon-quark coupling, which requires additional assumptions 
about the form of the nonlocal current, as discussed in Sec. [TT]. The resultant coupling in the 
present model and the means of calculating the electromagnetic decays of the vector mesons 
were discussed in Subsec. [Ill These photon- vector- meson couplings are defined by 

(o| jHp*> = -g^C 

(OlJ^Us) = -g^. (50) 

The empirical values, obtained from p — > e + e~ and to — > e + e~ j48|l , are g pl = 0.1177 GeV 2 
and = 0.0359 GeV 2 . The calculated values for these couplings, given in Table 5, are in 
reasonable agreement with the experimental ones. 

2 

Values for the dimensionless quantities gy = — - are also given in Table 5 and may be 

9v~{ 

compared with the results for g p7Tn in Table 4. Universal coupling of the p would predict that 



gp = g piviT (see, for example, |37|,|38[]). This is clearly violated in the model, although notably 
less so with parameter set A, where both of these couplings are closer to the empirical ones. 

It is interesting to compare the pj coupling with the corresponding coupling of the ai to 
the transverse axial current. A coupling strength g ai can be defined analogously to gy 7 , 

(0I4VJ = -g ai 5 ab ^ s . (5i) 

Evaluating this as described in Subsec. |III D| , we obtain values for g ai of 0.072 GeV 2 and 0.138 
GeV 2 with parameter sets A and B, respectively. No direct experimental measurement of g ai 
exists with which to make a comparison, but the quantity does appear in Weinberg's sum 



rules | 5q| . If one assumes complete vector and axial- vector meson dominance in Weinberg's 



first and second sum rules, then one gets the relations 

9 H - 4- = fl (52) 
m 2 p m 2 ai 

9(n = 9a x - (53) 

The results of the model for parameter set A are consistent with these vector- dominance 
versions of these sum rules, at the ~ 15% level. In contrast the results for set B clearly fail 
to satisfy these relations. 
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The pion form factor provides a further test of vector-meson dominance (VMD). This 
quantity receives contributions from the two kinds of diagram shown in Fig. 7. The first 
of these involves a triangular loop, and is often the only one considered in calculations of 
the form factor. For a timelike momentum q carried by the external current, the situation 
is similar to that discussed in Sec. [V] for the triangle diagrams in hadronic decays, with 
pseudo-threshold effects appearing at energies beyond twice the real part of the quark pole. 

The other kind of diagram in Fig. 7, which we shall refer to as a two-body diagram, 
is produced by the terms Gi(ij5 (g> 275) and ^2(7^75 <8> lu 7s) with type-II structure ([lOD in 
the nonlocal isovector current. The quark propagators involved are evaluated at p ± 
and so pseudo-threshold effects are not present. The contribution of these diagrams can 
be written as a sum over terms, each of which is the product of two loop integrals which 
are somewhat similar to those in Jpp, Jap or J\a- Since the path variable appears in form 
factors associated with both of the loops, the two loops are coupled together. The need 
for similar two-body contributions has also been noted in the context of models where the 
four-quark interaction is dependent only on the relative momentum of the qq pair |]IE|| . In 
a model of that type the analogue of the two-body diagram can be reduced to a single 
two-quark loop, where one of the the irqq vertex functions is modified by the photon. Such 
diagrams make no contribution to the pion charge, unlike the two-body diagrams in the 
present model. 

The absolute value of the form factor is displayed in Fig. 8, for a parameter set with 
m o(0) = 300 MeV. Below the p meson pole, very little variation with mo(0) is observed in the 
results when other parameter sets are used. Also shown on the figure are experimental data 
points, obtained from Refs. f5]]]5£| for the region of spacelike q 2 and Ref. |59j for timelike q 2 . 



Below the p pole, the model curve is in rather good agreement with the data, although its 
rise is a little shallower. This is confirmed by the mean-square pion radius, given in Table 5. 
The values are somewhat smaller than the experimental one [0] of 0.439 ± 0.008 (fm) 2 . 

The dashed curve in Fig. 8, labelled VMD, is plotted to test to what extent the model 
is consistent with the concept of vector meson dominance. Under the assumption of VMD 
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in the photon-pion coupling, the form factor is 



FM 2 ) = 1 " (54) 

p 



9 P T -m 



The ratio of g p7T7T to g p is underestimated by the model, and as a result the VMD form factor 
rises somewhat more slowly than the data. Nonetheless the VMD approximation to the 
model curve is not a bad one, particularly at low q 2 . 

An important check on our calculations, both analytic and numerical, is that the pion 
charge should be unity, ^(0) = 1. This result has been verified analytically in the chiral 



limit p2| . It is also satisfied by our numerical results, to the accuracy of our integration 
routines. About a third of the pion charge comes from the two-body diagram in Fig. 7, 
demonstrating the importance of including this contribution. 

The individual contributions to the form factor, below the p pole, are shown in Fig. 9. 
The curve labelled "bare" is the contribution from the triangle diagram with the local ^qq 
coupling (the first diagram in the representation of the dressed jqq vertex shown in Fig. 5) . 
The curve labelled "scalar" comes from the triangle diagram with a nonlocal coupling (as 
in the second diagram of Fig. 5) and a loop with a scalar insertion, given in Eq. ([$5|). This 
contribution is negligible over the range of q 2 considered. All other contributions are included 
in the curve labelled "vector" , corresponding to the part of Eq. ( fl3f) that is proportional to 
B(q 2 ). Since the transverse p-meson propagator is contained in the function B 2 (q 2 ) (fl4|), 
this contribution is the dominant one close to the p pole. The contribution of the two- 
body diagram is relatively small in the vicinity of the p. However it varies only slowly with 
momentum, and so it quite rapidly becomes important at large spacelike momenta. This is 
as expected since, for large momentum transfer to the pion, the pion vertex functions cut 
down the triangle-diagram contributions. 

It is interesting to note that, away from the pole, much of the variation with momentum 
is controlled by the "bare" piece of the form factor, which contributes ~ 77% of the mean- 
square charge radius (r 2 ). Although this term has no p pole, when added to the "vector" 
contribution, the result is close to that of the VMD approximation to the model. This 
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implies cancellation between the "bare" contribution and that of other states above the p in 
the "vector" piece, to leave the p pole as the dominant contribution. In contrast, although a 
similar mean-square radius has been obtained in an extended NJL model pOf , most of that 
value was accounted for by a diagram with an intermediate p meson, the bare photon vertex 
yielding just 32%. 

Above the p pole, the measured form factor is not well described by the model curve. 
A possible explanation is suggested by Fig. 10, which breaks the form factor down into its 
various contributions for q 2 > m 2 p . There is a substantial cancellation between the "bare" 
and "vector" contributions and so the result is likely to be very sensitive to fine details of 
the model in this regime. This is borne out by the strong dependence of the results in this 
region on the parameter set chosen. It is also responsible for the prominence of the rather 
strange structure just above the pseudo-threshold in Fig. 8. This cancellation between large 
amplitudes suggests that the results of the model should not be regarded as reliable in this 
region. 

The use of the conserved current constructed in Sec. |TI| implicitly ensures that the elec- 
tromagnetic Ward identities are satisfied by the model. Several examples of these have been 
already discussed above. Another important one is the amplitude for the decay n° — > 77, 
which is an example of an anomalous process. Such processes require a complete set of 
quark states and so present a problem for the usual NJL model, where the use of a reg- 
ulator means that high-energy quark states are discarded. In the model studied here, the 
low-energy theorem for tt° — > 77 is automatically satisfied provided that one includes both 
of the diagrams shown in Fig. 11. As well as the usual triangle diagram, there is a two-body 



diagram involving a dressed jqq vertex (|43| ) for one of the photons and a jqqqq vertex for 
the other. Details of the calculation are given in the Appendix. 

The finite current quark masses mean that the physical amplitude differs slightly from 
its value in the chiral limit. Numerical results for g^yy, defined from 

(7(?i)7(92)|7r ) = 2a f M g^e^^q" q%e^e? (55) 
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are given in Table 5. The deviations from |, the result in the chiral limit, are small and are 



consistent with those in the experimental value, = 0.503 ± 0.018 . 

The related process where one of the photons is off-shell, 77* — > tt°, provides a probe of 
the structure of the neutral pion. A corresponding form factor may be defined as 

As well as older data from the CELLO collaboration ||61|| , there is now new data from CLEO 



on this form factor in the spacelike region, which extends from — m\. These are shown 
in Fig. 12, along with the results of the model for a parameter set with m (0) = 300 MeV. 
The model results are not sensitive to the choice of parameter set and are dominated by the 
contribution from the triangle diagram with local photon couplings. They are consistent 
with the limited data available and comparable to the results obtained in other Bethe- 
Salpeter approaches |4T| , f42"fl . There are hopes that further experimental data will become 
available, since there are plans to measure the form factor with high precision at Jefferson 
Lab J53 . Also displayed on Fig. 12 is the VMD prediction for the form factor, given by 



f,~M 2 ) = 1 - — E — -r^> (57) 

#T77 yt^o, 9v q 2 -m 2 v 
using the model values for the couplings. (The couplings gvvp describing the decays V — > 
7T7, will be discussed shortly.) In this case, the VMD approximation to the model result is 
a little poorer than it was for the pion form factor. 

Finally, we consider decays of spin-1 mesons into 777 final states. Schematically, the 
calculation involves triangle and two-body diagrams, analogous to those of Fig. 11. Just as 
in the anomalous pion decay the piece of the nonlocal current which gives rise to the two-body 
diagrams for V — > 777 is the type-I term (||) with Dirac structure (^2(7^75 <8> 7^75)- These 
two-body diagrams do not prove to be numerically important in these decays, producing 
less than 1% of the total amplitudes. The results for these couplings, as defined by 



(7(91)71" (q 2 )\u) =ieg U7n e a p lu ,<ft<£ 



(7(gi)A<? 2 ) \p a ) = i5 ah e gpwj e a/3[iV q?q^e;, (58) 
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are given in Table 5, along with the corresponding decay widths. In the isospin symmetric 
case, there is no e ab3 component to the pir'-f matrix element, so that the decay widths for the 
charged and neutral p mesons are equal. The model results agree well with the experimental 



values PES 



T(co -> vr 7 ) = 717±43keV, 
r(p° -> 7r° 7 ) = 121 ±31keV, 

T(p ± -> vr ± 7 ) = 68 ± 8keV, (59) 

the difference between the measured charged and neutral p-decays not being considered 
statistically significant in view of the large error bars |)4] . 



Extending the uii'j amplitude to off-shell photon momenta, we compare our model results 
to the form factor as measured in Ref. p5| . The reaction 7t~p — > nu — > mr°p + p~ was used 
to probe the form factor in the range from 4m^ to (m^ — m^) 2 . Using a definition of 

the model results and experimental data are shown in Fig. 13. A parameter set where 
mo(0) = 300 MeV has again been used in plotting the model results, which are not sensitive 



to the set chosen. In common with other theoretical calculations [66] the model is in agree- 
ment with the data points at low q 2 but there is a discrepancy with the higher-energy values. 
It may be that there is some effect on this form factor due to the tail of the p' resonance 
f67| . Also, a potentially important missing ingredient in our calculation is ucj) mixing, since a 
calculation of this form factor within an SU(3) effective Lagrangian approach |JpJ has found 
a significant dependence on the mixing strength. Improved data would be needed to draw 
any firmer conclusions and there are hopes that the experimental situation will be indeed 
be clarified by forthcoming experiments at VEPP-2M or DA$NE [Q . 
In this case, the comparison with the VMD prediction, 

= 1 - ^ / 2 , (61) 
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is not completely straightforward, since the coupling g wpn cannot be calculated on-shell. 
Nonetheless by extrapolating to the soft-pion limit (zero pion four-momentum) we can get a 
reasonable estimate for the strength at the pole in Eq. (|6T|). The results for g W(m determined 
in this way are within 20% of the prediction of universal coupling, g^ pi , = g^^g p . The 
results for the form factor plotted in Fig. 13 again show good agreement between the full 
calculation and the VMD approximation. 

In the amplitude for the decay af — > 7r ± 7, the two-body diagram involves the terms in 



the nonlocal current Gi(ry5 ® 275) and (22(7" 75 ®7^75) which have a type-II structure (|10|). 
Since the a\ mass lies above the pseudo-threshold for parameter set B, in that case both the 
triangle and two-body diagrams require residue contributions. Gauge invariance imposes 
the following structure on the amplitude, 



^(qi)^ b (q 2 )\a a 1 )=te ab3 eg ai 



777 



2(q 2 ■ e ai )(q 2 ■ e;] 



(62) 



(m 2 ai - ml) 

The values calculated for <7 0l7ry , and the corresponding decay widths, are presented in Table 5. 
With parameter set A the dressed scalar piece of the nonlocal photon coupling and the two- 
body diagrams make relatively small contributions; for set B, these pieces are rather larger in 
magnitude but largely cancel with each other. The final results for all choices of parameters 



are much smaller than the experimental measurement [69] of 640 ± 246 keV 



VII. CONCLUSIONS 

We have investigated extended versions of the model proposed in Ref. J24| with nonlocal, 
four-fermion interactions, based on the separable form (|) suggested by instanton-liquid 



studies 1 15). The interaction form factors ensure the convergence of all quark-loop integrals 



and, for most of the acceptable range of parameters, also lead to quark confinement p4 |. 
The familiar problems of the NJL model are therefore eliminated whilst retaining much of 
the simplicity of that approach. 

Quark confinement is particularly important in applying such a model to heavier mesons 
than the pion, since it avoids a threshold for qq production occurring at an inconveniently 
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low energy, as in the NJL model. This feature of the nonlocal approach makes it especially 
well suited as the basis for the extended model that we study here, with interactions that 
bind vector and axial-vector mesons. 

The analytic structure of the quark propagator, with poles at complex momenta, means 
that a scheme has to be specified for continuing amplitudes to Minkowski space. We have 
used the scheme suggested by Lee and Wick |35| and Cutkosky et al. p6| , which leads to 



nonanalytic behaviour of the meson propagators above a pseudo-threshold. This behaviour 
ensures that no poles appear in the longitudinal vector channels. It also provides restrictions 
on both the admissible range of model parameters and the region of applicability of the 
model. 

In order to calculate the pion decay constant, which sets the basic scale for the model, as 
well as electromagnetic couplings, we need to construct conserved currents that are consistent 
with our nonlocal interaction. This has been done using the Noether-like method of Ref. 



24j. Ward identities related to the current conservation follow automatically and several 
have been analytically checked, including the Gell-Mann-Oakes-Renner relation and the 
low-energy theorem for 7r° — > 77 in the chiral limit. The latter involves the axial anomaly, 
which has formed a long-standing problem for the usual NJL model with its cut-offs on the 
quark propagators. 

Working at tree level in terms of mesons (leading order in 1/N C ) and to all orders in 
momentum, we have calculated the meson masses and various strong decay rates. When the 
p mass is used to fix the strength of the relevant four-quark vertex, we find that the p-meson 
width is reasonably well described. In particular the model provides an improvement on the 
underestimated value obtained in the extended NJL model ||50|| . As in the NJL model, a 



light a is found, with a mass similar to that found in a number of analyses of tttt scattering 
The calculated mass of the a\ meson is somewhat lighter than the observed value. By 
cutting down the available phase-space, this also means that the model gives too small a 
width for the decay a\ — > pn. 

Various electromagnetic quantities have also been calculated in this model. The decay 
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rates found in the model are in reasonably good agreement with the observed ones, except 
for the case of a± — > Try. The electromagnetic form factor of the pion agrees well with 
the data, at least below the p pole. The form factors for 77* — > n and uj — > ny* are also 
in agreement with the, admittedly rather limited, data currently available. We have also 
compared our results for electromagnetic form factors with vector- dominance formulae using 
on-shell couplings as calculated in the model. Although diagrams where the photon couples 
directly to the quarks make a significant contribution, they combine with diagrams that 
involve an intermediate vector meson to give results that are very close to those of VMD, 
at least at low momenta. The model is thus able to provide some support for the idea of 
vector dominance in photon-meson interactions. 

Overall the level of agreement with observed meson properties is satisfactory, given the 
simple nature of the model and the fact that important effects have not been included at the 
current level of approximation. It is desirable to extend the current work to include meson- 
loop effects (l/N c corrections) since these could be significant. In particular, they may be 
important for the scalar isoscalar channel, where the a meson is very strongly coupled to 
two-pion states and where there are some indications that our leading-order calculations do 
not agree with phenomenological analyses. 
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APPENDIX A: ir° -► 77 

We outline here the various contributions to the n° — > 77 amplitude in this model. We 
work in the chiral limit, where the axial anomaly requires that g nll = \. The relevant 
diagrams are shown in Fig. 11. Consider first the simpler case where vector meson degrees 
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of freedom are not present in the model, setting G 2 = G3 = 0. The triangle diagrams where 
both of the photons couple to the local current, 7 M , contribute 



^e a ^q p 2 eTel l 'e 2 g nqq J J _ J I -2m (p 2 ) + 4m' (p 



2 



(2tt) 4 [p 2 - m 2 (p 2 )Y V q 

+ 4ra' o{p ^^-"f ) , (A!) 

(91 - 12) 2 J 

where the prime denotes a derivative with respect to p 2 . The contribution of diagrams where 
one of the photons couples to the nonlocal piece of the current, with a scalar loop insertion 

Diagrams with a nonlocal coupling at both photon vertices do not contribute to the ampli- 
tude because the resulting Dirac trace vanishes. Converting the integral in the sum of ( |A1[ ) 
and ( |A2| ) into Euclidean space, and changing variable to 

t = (A3) 
Pe 



leads to 



7T m (0) JO (1 + £) d 

The low-energy theorem now follows by invoking the analogue of the Goldberger-Treiman 



P *ix *v a EM 9-Kqq f°° dt 



relation 70 in this model, 



9-Kqq — (^5) 
In 

which was derived in Ref. |24]] by considering f n in the chiral limit .0 

In the extended model with vector mesons, g nqq is affected by the pseudoscalar-axial 
mixing induced by the G 2 coupling ( |2"8| , |29"|) . However, the form of f n is also modified, as 
discussed in Subsec. [Ill D| , and hence Eq. ( |A5| ) remains valid in the chiral limit ||32j| . In 



'In the notation of Ref. [ 24 1 g nqq is Z£ ■ 
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addition the photons are on-shell and so we need only consider the bare 'jqq vertex and the 
nonlocal piece fl35|) with a scalar one-quark loop, as also discussed in Subsec. [Ill D| . The sum 
of the contributions (|A1| ) and ( |A2| ) therefore yields g wll = |, just as in the model without 



vector mesons. 



This feature is in agreement with the analysis of Ref. [|7TJ, where it was shown that for 
a quark propagator without wave function renormalization, then the anomaly is saturated 
by taking only the leading part of the pion Bethe-Salpeter amplitude and the dressed jqq 
structures subject to the Ward identity for that vertex. The statement is non-trivial because 
terms in the pion amplitude that are linear in momentum can contribute to the decay 
amplitude, even in the chiral limit. For instance, the ^75 term which appears in V^q) for 
this model gives rise to additional triangle-diagram contributions. From such diagrams with 
two local photon vertices, one finds 



m n J (2vr) 4 [p 2 - ml(j> 2 )Y \ * ^' q 2 

+ ^: {qi - q f ) . <A6) 

Since g nqq is of 0(m n ) in the chiral expansion, this expression is of 0(1). Similar diagrams 
with one local and one nonlocal vertex give 



p *a< *u 29n qq f d A p f 2 (p)8mo(p 2 )m' (p 2 ) ( 2 (p ■ qf [p ■ (gi - q 2 )} 2 ' 



Me^q^e, e 2 e ^ J \p*-ml{p 2 )Y V q 2 (ft - g 2 ) 2 



(A7) 



(The diagrams with two nonlocal vertices again have a vanishing Dirac trace.) 

Since the sum of ( |A~6| ) and ( |A7|) is non-zero, there must be some contribution that cancels 



them in the full amplitude for the anomalous process. The relevant piece arises from the 
two-body diagram, which is displayed in Fig. 11. Terms in the nonlocal vector current of 
the form G^T^Ts ® 7^75) with a type- 1 @ structure are responsible for the jqqqq vertex. 
(Note that for simplicity we have suppressed isospin factors.) This diagram factorizes into 
two separate loop integrals, the loop between the two photons giving rise to the e Q( g Mi , factor. 
The other loop is simply a linear combination of the familiar integrals JAp(tTi 2 ) and J^ A (m 2 ) 
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from Subsec. 1 1 1B|. This combination can be simplified by recalling the definitions of the 
pion-quark coupling constants in (p8|). In the chiral limit, the contribution of this diagram 
can be written as 

J (2vr) 4 [p 2 - m 2 (p 2 )\ 2 \ q 2 (gi - g 2 ) 2 / 

Converting to Euclidean space and integrating by parts, this expression can be shown to 



cancel exactly with the sum of flAjp and (|A7|) , demonstrating that the low-energy theorem 
for 7T° — > 77 holds in the extended model. 
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APPENDIX: TABLES 



Parameter 


Set A 


Set B 


Parameter 


Set A 


Set B 


m (0)(MeV) 


280 


320 


m(0)(MeV) 


326 


370 


m c (MeV) 


8.4 


11.0 


A(MeV) 


995 


846 


Gi(GeV- 2 ) 


37.1 


57.6 


G 2 (GeV- 2 ) 


-5.70 


-6.53 


G 3 (GeV- 2 ) 


-5.20 


-5.86 


G 4 (GeV- 2 ) 


-0.80 


-4.14 


G 5 (GeV- 2 ) 


2.57 


4.76 









Table 1. Values of the model parameters for sets A and B, fitted as discussed in the text. 
Also shown is the dynamical quark mass m(0). 



Set A 


Set B 


±496 ± 130i 
±1168 ± 790i 
±1488 ± 1155i 
±1742 ± 1436i 


±404 ± 257i 
±962 ± 702i 
±1242 ± 1005i 
±1463 ± 1240i 



Table 2. Positions of the lowest four sets of poles in the quark propagator. 
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Set A 


Set B 


Particle 


Mass 


9iqq 


9iqq 


Mass 


9iqq 


9iqq 


7T 


Fit 


3.44 




0.0739 


Fit 


3.91 


0.0715 


a 


443.2 


3.51 




- 


465.8 


4.06 


- 


P 


Fit 


1.12 




- 


Fit 


1.11 


- 


d\ 


946.8 


1.13 






1061.5 


2.27 


- 


to 


Fit 


1.07 






Fit 


1.05 




h 


Fit 


0.89 






Fit 


2.51 




a 


Fit 


0.75 






Fit 


1.71 




rf 


874.9 


0.83 




0.190 


899.4 


2.36 


1.448 



Table 3. Meson masses (in MeV) and couplings of the mesons to quarks. The couplings gi qq 
and g iqq are defined by the equations fl2"?D and (|28|). 





Set A 


Set B 


Coupling 


Value 


Width(MeV) 


Value 


Width(MeV) 


g a ^(MeV) 


1438 


108.0 


1625 


135.1 


9p7V7T 


5.52 


126.0 


5.26 


114.0 




10.65 


74.0 


11.77 


116.4 


g aifm (MeV) 


2174 


44.0 


4604 


376.2 




18.19 




10.87 




R 


-0.048 




-0.087 





Table 4. The meson couplings, as defined in equation (^), together with the corresponding 
partial widths. R is the ratio of the d- to s-wave amplitudes defined in ([491). 
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Set A 


Set B 




Set A 


Set B 


Quantity 


Value 


Value 


Quantity 


Value 


Value 


S P7 (GeV 2 ) 


0.0889 


0.0773 


9 P 


6.67 


7.67 


WGeV 2 ) 


0.0308 


0.0265 


9u> 


19.92 


23.12 


<r 2 )(fm 2 ) 


0.346 


0.344 


- 


- 


- 




0.505 


0.501 


- 


- 


- 


SW 7 (GeV _1 ) 


-2.29 


-2.25 


T(lo^ 7r 7 )(keV) 


692 


669 


^(GeV- 1 ) 


-0.755 


-0.707 


T(p -> 7r 7 )(keV) 


71.6 


62.7 




140.2 


201.5 


r(ai -> tt 7 ) (keV) 


24.7 


45.7 



Table 5. Electromagnetic couplings of mesons. 
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APPENDIX: FIGURES 




Fig. 1. Schwinger-Dyson equation for the quark propagator in the ladder approximation. 




Fig. 2. Bethe-Salpeter equation for qq scattering in the ladder approximation. 





□t 


P4 














\ • 


• 1 





Fig. 3. Deformed integration contour in plane, beyond the pinch point. The open boxes 
denote the poles of the p_ propagator and the filled circles those of the p + propagator. The 
arrows indicate the directions in which these poles move as go increases. 
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Fig. 4. Coupling of particle to an external current. V denotes the vertex function fl25p. 




Fig. 5. The dressed 'yqq vertex. T denotes the qq scattering matrix in either the transverse 
or longitudinal vector channel. 
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Fig. 6. Denominator of the propagator in the sigma channel, 1 — G\Jss, and the pion 
determinant, (^), as functions of the timelike meson momentum. The denominators of the 
p and <2i propagators, 1 — G2 Jvvaa^j are & l so displayed, scaled by a factor of 10. 
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Fig. 7. Spacelike pion form factor. There is another similar triangle diagram where the 
photon couples to the anti-quark. 
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Fig. 8. Absolute value of the pion form factor, plotted against q 2 in GeV 2 . The solid 
line is the model result; the dashed line is the VMD approximation to it. At the level of 
approximation used the form factor has a pole at the p-meson mass. The data points are 
from Refs. 0-0]. 



51 




52 



6 







-2 



-8 



-10 



- 




- 


Bare 




Two- body 




Scalar 


- 


/ Fu form factor 


: , / 


/ Vector 



0.6 0.65 0.7 0.75 0.8 0.85 0.9 0.95 1 



2 

q 



Fig. 10. Various contributions to the timelike pion form factor, above the p pole, plotted 
against q 2 in GeV 2 . The different contributions are defined in the text. 
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Fig. 11. 7T° — > 77. There are also similar diagrams with a <-> 6. 
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Fig. 12. The 7177* (q 2 ) transition form factor, plotted against q 2 in GeV 2 . The solid line is 
the model result; the dashed line is the VMD approximation to it. The open data points 
are from Ref. pl |, the solid ones from Ref. 
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Fig. 13. The square of the c<j7r7*(g 2 ) form factor, plotted on a logarithmic scale against q 2 
in GeV 2 . The solid line is the model result; the dashed line is the VMD approximation to 
it. The data points are from Ref. 
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